Thermal entanglement in the anisotropic Heisenberg XXZ model with the Dzyaloshinskii-Moriya 

interaction 
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The entanglement is investigated in a two-qubit Heisenberg XXZ system with Dzyaloshinskii-Moriya(DM) 
interaction. It is shown that the entanglement can be efficiently controlled by the DM interaction parameter and 
coupling coefficient Jz . Dx (the x-component parameter of the DM interaction) has a more remarkable influence 
on the entanglement and the critical temperature than Dz(lhe z-component parameter of the DM interaction). 
Thus, by changing the direction of DM interaction, we can get a more efficient control parameter to increase the 
entanglement and the critical temperature. 
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I. INTRODUCTION 

Entanglement has been extensively studied in recent years 
because it has the fascinating nonclassical nature of quantum 
mechanics, and it plays a key role in quantum information 
processing Jll |31- The quantum entanglement in solid state 
systems such as spin chains is an important emerging field p, 
BSIal^ spin chains are natural candidates for the realization 
of entanglement, and spin has been researched in many other 
systems, such as superconductor 13, S], quantum dots 1191 UOl 
Qj, and trapped ion ifllfisll. 

In order to characterize qualitatively and quantitatively the 
entanglement properties of condensed matter systems and ap- 
ply them in quantum information, the thermal entanglement 
qualities in Heisenberg model have been extensively studied 
11141 US. [la, UTl [l8ll . and many schemes of teleportation via 
thermal entangled state have been reported 11191 l20l un. l22l 
I23II . In condensed matter systems, the Heisenberg chains have 
also been used to construct a quantum computer ll24ll . perform 
quantum computation ll25ll26ll27ll . and so on. 

In those studies the spin-spin interaction was considered, 
but the spin-orbit coupling was rarely considered. Espe- 
cially the influences of x-component DM interaction param- 
eter on the entanglement and the critical temperature have 
never been reported. In this paper we investigate the influence 
of Dzyaloshinskii-Moriya interaction parameter (arising from 
the spin-orbit coupling) and coupling coefficient J^ on the en- 
tanglement of a two-qubit anisotropic Heisenberg XXZ spin 
chain. We show that the DM interaction parameter and the 
coupling coefficient J^ are both efficient control parameters 
of entanglement, increasing them can enhance the entangle- 
ment or slow down the decreace of the entanglement. In addi- 
tion, by analyzing we know that different component param- 
eters of DM interaction have different influences on the en- 
tanglement and the critical temperature Tc, the parameter D^ 
(x-component parameter of the DM interaction) has a more 



remarkable influences than the parameter D^ (z-component 
parameter of the DM interaction). So more efficient control 
parameter can be obtained by changing the DM interaction 
direction. 

Our paper is organized as follows. In Section II, we intro- 
duce the Hamiltonian of the two-qubit anisotropic Heisenberg 
XXZ chain with the z-component parameter of the DM inter- 
action, we calculate the concurrence of this system, analyze 
the influence of parameters on the entanglement in ground 
state and thermal state. In Section III, we similarly analyze 
the model of the two-qubit Heisenberg XXZ chain with the 
x-component parameter of the DM interaction. Then we com- 
pare the influences of the two component parameters of the 
DM interaction on the entanglement in Section IV. Finally, in 
Section V a discussion concludes the paper. 



II. XXZ HEISENBERG MODEL WITH DM 
INTERACTION PARAMETER D^, 

The Hamiltonian H for a two-qubit anisotropic Heisenberg 
XXZ chain with DM interaction parameter D^ is 



H = Jcrfcr^ 



Jcflal + J^ala'^+D^ialcjl-cjlcj^), (1) 



where J and J^ are the real coupling coefficients, D^ is the z- 
component parameter of DM interaction, (T^{i = x, y, z) are 
Pauli matrices. The coupling constants ,/ > and J^ > 
correspond to the antiferromagnetic case, J < and J^ < 
correspond to the ferromagnetic case. This model is reduced 
to isotropic XX model when J^ = and isotropic XXX model 
when Jz = J- Parameters J, J^ and D^ are dimensionless. 

In the standard basis {|00), |01), |10), |11)}, the Hamilto- 
nian (1) can be expressed as 



H = 
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By a straightforward calculation we can obtain H eigenstates: 
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\^i) - |oo), 



(3a) 



1*2) = 111) 



|$3) = ^(e^'|01) + |10)), 



|<I>4) = -^(e^^|01)-|10)), 

with corresponding eigenvalues: 

El ~ Jz, 

E2 ~ Jz, 

E3 = -Jz + 2w, 



(3b) 
(3c) 

(3d) 

(4a) 
(4b) 
(4c) 




FIG. 1 : The concurrence is plotted versus T and Dz where the cou- 
pling constants J — 1 and Jz = 0.2. 



E4 = —Jz — 2w, 



(4d) 



where w = yO^TO^, 9 = arctan(-^). 

The state of a spin chain system at thermal equilibrium is 
p{T) = °^P(-^g) , where Z = tr[exp(-/3H)] is the parti- 
tion function of the system, H is the system Hamiltonian and 
f3 — K^ T ' "^^^^ T is temperature and Kb is the Boltzmann 
costant which we take equal to 1 for simplicity. Here p{T) 
represents a thermal state, so the entanglement in the thermal 
state is called thermal entanglement 12811 . In the above stan- 
dard basis, the state of this system at thermal equilibrium can 
be expressed as; 



piT) ^ - 
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where 



u 



iPw\MJ,-2w) 



')e 



i(l + e4/3-)e'3(^^-2-), y = 1(1 
and Z = 2e~'^'^- [1 + e^^'^- cosh(2/3w)]. 



In what follows, we consider the concurrence to quantify 
the amount of entanglement of the above two-qubit system 
state /3(T). The concurrence 1129113011 is defined as C{p{T)) — 
max[2max(Ai) — S^AijO], where \i{i — 1,2,3,4) are the 
square roots of the eigenvalues of the matrix R = pSp* S, in 
which S = a^ ^ af , p is the density matrix of the Eq. (5), 
and the asterisk denote the complex conjugate. After some 
straightforward calculation, we get 



1 
Z~ 



Al=4e-''^^ 



A3 = 4e''('^^-2^), 



A4 - ^e 



(6a) 



(6b) 



(6c) 



(6d) 



thus the corresponding concurrence can be expressed as: 



CipiT)) 



^{- 



^PJ. 



^2Pw_ 



-2/3J,| 



-2/3 1«_ 



z 



-2^-'^),0}. 

(7) 

The concurrence is invariant under the substitutions J ^ —J 
and Dz -^ —Dz, so we can restrict J > and Dz > 
without loss of generality. The concurrence ranges from 
to 1, C{p{T)) = and C{p{T)) = 1 indicate the vanish- 
ing entanglement and the maximal entanglement respectively. 
We can see from Eq. (7) that the entanglement C(p{T}) = 



(e2/3^ _ e- 



-2/}./. _ „-2/3u) _ „-2/3J, 



) if J^ > — w, and 
C{p{T)) = if Jz < — u'. Here we analyze the Jz > —w 
case. 

When r = 0, the system is in its ground state. It is easy to 
find that the ground-state energy is equal to 



E4 = -Jz -~'^w,ifJz > 



El ^ E2^ Jz.ifJz < -w. 



(8a) 



(8b) 



Thus, the ground state is the disentangled state l^i) or |$2) 
when Jz < —w, and the ground state is the entangled state 
I $4) when Jz > —w. The entanglement of the ground state 
|$4) is the maximal entanglement with C(|$4)) = 1. 

As the temperature increases, thermal fluctuation will be 
introduced into the system, thus the entanglement will be 
changed due to the mix of the ground states and the excited 
states. When the temperature is higher than a critical tem- 
perature the entanglement is zero. Quantum phase transition 
happens at the critical temperature Tc- From Eq. (7), we ob- 
tain the following critical temperature equation 



smh(— ) 



(9) 



To see the change of the entanglement in detail, we analyze 
the concurrence of Eq. (7). By fixing some parameters we 
can know the roles of the other parameters and the variation 
of the entanglement. In Fig. 1, the thermal entanglement is 




FIG. 2: The concurrence is plotted as a function of tiie temperature 
T for different D^, here J = 1 and J^ = 0.2. 




glement or increase the critical temperature to slow down the 
decreace of the entanglement. 



III. XXZ HEISENBERG MODEL WITH DM 
INTERACTION PARAMETER D^ 

Here we consider the case of the two-qubit anisotropic 
Heisenberg XXZ chain with DM interaction parameter D^- 
The Hamiltonian is 

h' ^Jaf(J^ + Jcr\a« + Jz(rt<^2+Dx{<jl<J2-<^l<^2)^ (10) 

where D^ is the x-component parameter of DM interaction, 
J, Jz and (j'^{i = x, y, z) are the same as those in Section II. 
Parameters D^, J and J^ are dimensionless. 

In the standard basis {|00), |01), |10), |11)}, the Hamilto- 
nian (10) can be rewritten as 



H = 



Jz iD., -iD, 

—iDx —Jz 2 J iDx 

iDx 2 J —Jz —iDx 

—iDx iDx Jz 



(11) 



After a straightforward calculation we obtain H eigenstates: 

1 



|vl/i)=.— (|00) + |11)), 



|VI/2) = ^(|01) + |10)), 



(12a) 



(12b) 



FIG. 3: The concurrence versus T for different J^ in the system with 
z-component parameter D~, here J = 1 and D^ = 1. 



plotted versus T and Dz where the coupling constants J = 1 
and Jz = 0.2. From the figure, it's obvious that the increased 
temperature T can decrease the entanglement. The reason is 
the mixing of the maximally entangled state with other states. 
In addition, it is easy to see that the entanglement will increase 
as the DM interaction parameter Dz increases. Fig. 2 demon- 
strates the concurrence versus temperature for different DM 
coupling parameter Dz when J = 1 and Jz = 0.2. It shows 
that the concurrence will decrease with increasing tempera- 
ture T and increase with increasing Dz for a certain temper- 
ature. The critical temperature Tc determined by Eq. (9) is 
dependent on Dz, increasing Dz can increase the critical tem- 
perature above which the entanglement vanishes. Similarly, 
Fig. 3 shows the concurrence versus temperature for different 
anisotropic coupling parameter Jz when J ~ 1 and Dz — I- 
We can see that the entanglement decrease with the increase 
of temperature, and by increasing Jz, the critical temperature 
is increased and the entanglement is enhanced for a certain 
temperature. 

So the DM interaction parameter Dz and anisotropic cou- 
pling coefficient Jz are both efficient control parameters of 
entanglement, by increasing them, we can enhance the entan- 



|*3> = ^(-isin0|OO)+cos(/)|Ol)-cos(/)|lO)+isin(/)|ll)), 
v2 

(12c) 

1*4) = ^(-isin.^|00)+cos<^|01)-cos(/3|10)+isin(^|ll)), 
v2 

(12d) 

with corresponding eigenvalues: 



El = Jz 



E2 — 2J — Jz, 



E^ = -J + w , 



E, 



-J — w , 



2D^ 



J+J^+w 



(13a) 



(13b) 



(13c) 



(13d) 
-), and 



where = arctan( , , 'i / ), (p ~ arctan( 

w' = ^{J + JzY + ^Dl 

For convenience of analysis, we assume Jz ^ J in this 
section (for the case Jz > J we can get some similar results). 
Here, it is easy to see that the system ground-state energy is 
E^ = —J — w , thus the corresponding ground state {'^'4) is 
the maximally entangled state with C(|^4)) = 1. In fact the 



four eigenstates in Eq. (12) are all the maximally entangled 
states, this phenomenon indicates that the ground states have 
more entanglement in this system than the system of Eq. (1). 
At thermal equilibrium the density matrix of this two-qubit 
spin chain system has the following form: 



P{T) = 



exp[~[3H ) 



2Z' 




,(14) 



where Z' = ^e'l^^ cos\i[l3{J - J^)] + 2e'^'^cosh(/3u;') is 
the partition function of the system, H is the system Hamil- 
tonian and /3 = ttT ^^'^^ ^^^ Boltzmann costant Kb = 



v± 



.I3(J.~ 



-2J) ± (, 



MJ-w ) 



cos^ (f) + e'^''^+'" ) cos^ if), and 



o/3(./-«' ) , 



^PiJ+w 



^ ^ le^^" ""'" ' sm(^C0S(^ + le^"'-" ''"" ' smipcosip. 

In what follows, we calculate the square roots of the eigen- 
values of the matrix R = p Sp *S, where p * is the complex 
conjugate of p and S = cti (S) o"! ■ The square roots of the 
eigenvalues of the matrix R are; 



Ai 



1 

— T^ 

z 



,/3(J,-2J) 



a; = ^re-^J' 



(15a) 



(15b) 



A3 = —e/''[cosh{Pw ) + \/cosh^{(3w') - 1], (15c) 



\ = — e^-^[cosh(/3w') 

Zj 



cosh2(/3w')-l]. (15d) 



According to the methods in II29L 13011 . when Jz ^ J we obtain 
the corresponding concurrence 



C(p'(T)) = max{47[|e^''(cosh(/3u'') + Vcosh^ (/3u;' ) - 1) 

Zj * 

- e-l^^--\ ~ e^^coshiPw) - ^J cosh^ {Pw' ) ~ I) 

- e'^(^^-2-^)],0}. (16) 

To analyze the role of parameters and the variation of the 
entanglement, we restrict D^ > 0, Jz > and J > 0. In 
Fig. 4, the concurrence is plotted versus T and D^ when the 
coupling constants J = 1 and J^ = 0.2. It is evident that 
increasing temperature will decrease the entanglement, and 
increasing Dx will enhance the entanglement and increase the 
critical temperature T^ which is determined by Eq. (16). Fig. 
5 demonstrates that the concurrence versus T for different Jz 
with x-component parameter D^ = 1 and J = 1, it is easy 
to find that increasing Jr can increase the critical temperature 
and enhance the entanglement for a certain temperature. So 
Dx and J^ are both efficient control parameters of entangle- 
ment, too. 




FIG. 4: The concurrence is plotted versus T and D^ where the cou- 
phng constants J = 1 and Jz = 0.2. 
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FIG. 5: The concurrence versus T for different Jz in the system with 
x-component parameter Dx, here J = 1 and D^ = 1. 



IV. THE COMPARISON BETWEEN THE TWO DM 
INTERACTION COMPONENT PARAMETERS 

From Section II and Section III, we know that the x- 
component parameter D^ and the z-component parameter D^ 
of DM interaction have the similar qualities. They are both 
efficient control parameters of entanglement, increasing them 
can enhance the entanglement or increase the critical temper- 
ature to slow down the decrease of entanglement. 

In this section, we mainly analyze the differences between 
x-component parameter and z-component parameter. We note 
that there is less disentanglement region in Fig. 4 than Fig. 1, 
where J — \ and Jz = 0.2, and the temperature and the spin- 
orbit coupling parameter have the same range. We also see 
that increasing x-component parameter D^ can make the en- 
tanglement increase more rapidly, for example, when T = 6 
the concurrence increases more rapidly in Fig. 4 than Fig. 1. 
These phenomena denote that x-component parameter D^ has 
a more remarkable influence than the z-component parameter 
Dz- In Fig. 6, the concurrence is plotted as a function of 
the temperature T for Dz — 2 and D^ = 2 with J — 1 




gated. The DM interaction parameter and coupling coefficient 
Jz are efficient control parameters of the entanglement. By in- 
creasing the parameters, we can enhance the entanglement or 
increase the critical temperature to slow down the decreace of 
the entanglement. In addition, we have also investigated the 
differences between x-component DM interaction parameter 
Dx and z-component parameter D^- Entanglement can be in- 
creased more rapidly by increasing D^ than D^- When D^ 
and Dz have the same value, D^ has a higher critical tem- 
perature than Dz- Thus, by changing the direction of DM 
interaction, we can obtain a more efficient control parameter 
to increase the entanglement and the critical temperature. 



FIG. 6: The concurrence is plotted as a function of the temperature 
T for _Dz = 2 and D^ = 2, here J = 1 and J^ = 0.2. 



and Jz — 0.2. It is easy to see that for the same Dx and 
Dz, x-component parameter D^. has a higher critical temper- 
ature and more entanglement for a certain temperature than z- 
component parameter Dz- We show distinctly the differences 
between different DM interaction parameters in Fig. 6. So 
using the DM interaction parameters in different directions, 
we can increase the entanglement and the critical temperature 
with different efficiencies. 



V. DISCUSSION 

The entanglement of a two-qubit Heisenberg XXZ system 
with different DM interaction parameters have been investi- 
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